2205.02501v2 [astro-ph.EP] 11 May 2022 


. 
. 


arXiv 


DRAFT VERSION Mav 12, 2022 
Typeset using IATEX twocolumn style in AASTeX63 


The Discovery of a Planetary Companion Interior to Hot Jupiter WASP-132 b 


BENJAMIN J. Honp,^?? KNICOLE D. CoróN,?? TRAVIS A. BERGER,”?'* VESELIN Kosrov,?? 

MICHELE L. SitvERSTEIN,??:* KEIVAN G. StTassun,* JACK J. LISSAUER,? KAREN A. COLLINS, RICHARD P. SCHWARZ,’ 
RAMOTHOLO SEFAKO,® CARL ZIEGLER,’ CESAR BRICENO,!? NicHoLAs Law, ANDREW W. MaNN,!! GEORGE R. RICKER,” 
Davip W. LATHAM, S. SEAGER,? 4 JosuuA N. Winn,’ Jon M. JENKINS,? LUKE G. Bouma,’© BEN Fark," 

GUILLERMO TORRES, JoskPH D. TWICKEN,* 5 AND ANDREW VANDERBURG!’ 


! Department of Astronomy, University of Maryland, College Park, MD 20742, USA 
2 NASA Goddard Space Flight Center, 8800 Greenbelt Road, Greenbelt, MD 20771, USA 
3 GSFC Sellers Exoplanet Environments Collaboration 
^ Department of Physics and Astronomy, Vanderbilt University, Nashville, TN 37235, USA 
5NASA Ames Research Center, Moffett Field, CA 94035, USA 
6 Center for Astrophysics | Harvard & Smithsonian, 60 Garden Street, Cambridge, MA 02138, USA 
T Patashnick Voorheesville Observatory, Voorheesville, NY 12186, USA 
8 South, African Astronomical Observatory, P.O. Box 9, Observatory, Cape Town 7935, South Africa 
9 Department of Physics, Engineering and Astronomy, Stephen F. Austin State University, 1936 North St, Nacogdoches, TX 75962, USA 
10 Cerro Tololo Inter-American Observatory/NSF's NOIRLab, Casilla 603, La Serena, Chile 
11 Department of Physics and Astronomy, The University of North Carolina at Chapel Hill, Chapel Hill, NC 27599-3255, USA 


1? Department of Physics and Kavli Institute for Astrophysics and Space Research, Massachusetts Institute of Technology, Cambridge, MA 
02139, USA 


13 Department of Earth, Atmospheric and Planetary Sciences, Massachusetts Institute of Technology, Cambridge, MA 02139, USA 
14 Department of Aeronautics and Astronautics, MIT, 77 Massachusetts Avenue, Cambridge, MA 02139, USA 
15 Department of Astrophysical Sciences, Princeton University, 4 Ivy Lane, Princeton, NJ 08544, USA 
16 Cahill Center for Astrophysics, California Institute of Technology, Pasadena, CA 91125, USA 
17 Space Telescope Science Institute, 3700 San Martin Drive, Baltimore, MD, 21218, USA 
18 SETI Institute, Mountain View, CA 94043, USA 


ABSTRACT 


Hot Jupiters are generally observed to lack close planetary companions, a trend that has been 
interpreted as evidence for high-eccentricity migration. We present the discovery and validation of 
WASP-132 c (TOI-822.02), a 1.85 + 0.10 Rẹ planet on a 1.01 day orbit interior to the hot Jupiter 
WASP-132 b. Transiting Exoplanet Survey Satellite (TESS) and ground-based follow-up observations, 
in conjunction with vetting and validation analysis, enable us to rule out common astrophysical false 
positives and validate the observed transit signal produced by WASP-132 c as a planet. Running the 
validation tools vespa and TRICERATOPS on this signal yield false positive probabilities of 9.02 x 107° 
and 0.0107, respectively. Analysis of archival CORALIE radial velocity data leads to a 3c upper limit 
of 28.23 ms~! on the amplitude of any 1.01-day signal, corresponding to a 3e upper mass limit of 37.35 
Ma. Dynamical simulations reveal that the system is stable within the 3e uncertainties on planetary 
and orbital parameters for timescales of ~100 Myr. The existence of a planetary companion near 
the hot Jupiter WASP-132 b makes the giant planet's formation and evolution via high-eccentricity 
migration highly unlikely. Being one of just a handful of nearby planetary companions to hot Jupiters, 
WASP-132 c carries with it significant implications for the formation of the system and hot Jupiters 
as a population. 

1. INTRODUCTION Ever since the Nobel Prize-winning discovery of the 

first exoplanet around a Sun-like star (Mayor & Queloz 

1995), hot Jupiters have represented one of the great- 

Corresponding author: Benjamin J. Hord est enigmas of exoplanet science. With radii Rp > 8 Re 
benhordGastro.umd.edu and orbital periods P « 10 days (Winn et al. 2010; Wang 
et al. 2015; Garhart et al. 2020), hot Jupiters represent 
a class of planets with no analogue in our solar system. 
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Traditional theories on planet formation are insufficient 
to explain the existence of giant gaseous planets so close 
to a host star (Lin et al. 1996). Therefore, new forma- 
tion scenarios have been put forth (e.g. disk migration, 
planet-planet scattering, secular migration) to explain 
the existence of hot Jupiters, most of which involve an 
inward migration after initially forming beyond the ice 
line (e.g., Lin et al. 1996; Rasio & Ford 1996). However, 
none of these mechanisms alone can satisfy all observ- 
able constraints, leaving the primary pathways of hot 
Jupiter formation largely still uncertain (Dawson et al. 
2014; Dawson & Johnson 2018). 

One clue that may help distinguish between different 
formation pathways is that hot Jupiters are predomi- 
nantly “lonely”, meaning they are the only planet in 
their system within a factor of 2 or 3 in orbital dis- 
tance (e.g. Steffen et al. 2012; Knutson et al. 2014; Endl 
et al. 2014; Hord et al. 2021), although they may have 
more distant companions, particularly giant companions 
(Schlaufman & Winn 2016). This lack of nearby com- 
panions is expected from a high-eccentricity migration 
formation pathway - a scenario in which the hot Jupiter 
migrates inwards from beyond the ice line via some form 
of gravitational perturbations that put it on an eccen- 
tric orbit that eventually circularizes much closer to the 
host star (Rasio & Ford 1996). This high-eccentricity 
migration results in the scattering and possible ejection 
of other planets in the system as the hot Jupiter's eccen- 
tric orbit sweeps through the inner parts of the stellar 
system (Mustill et al. 2015). 

Of the ~500 hot Jupiters currently confirmed, only 
three have proven to be exceptions to this lonely trend so 
far. The systems WASP-47 (Becker et al. 2015), Kepler- 
730 (Canas et al. 2019), and TOI-1130 (Huang et al. 
2020) all host a hot Jupiter with at least one nearby 
companion planet, making the high-eccentricity migra- 
tion scenario for these planetary systems unlikely!. It 
is more likely that these planetary systems formed via 
disk migration, where the protoplanets migrate inwards 
all together within the disk, potentially preserving plan- 
ets near the hot Jupiter (Lin et al. 1996; Lee & Peale 
2002; Raymond et al. 2006). These systems serve as 
rare opportunities to dynamically constrain the forma- 
tion of hot Jupiters and also potentially serve as a 


l We note WASP-148b is a hot Jupiter with an outer, massive 
companion at just within 3 times the orbital distance of the hot 
Jupiter. We choose not to include this system in our discussion of 
hot Jupiters with nearby companions because of its very different 
architecture from the other systems (which have smaller, closer 
companions) and because the mass and orbital distance of the 
companion planet puts it on the borderline of the approximate 
definition of hot Jupiters with close planetary companions. 


bridge to the slightly cooler population of warm Jupiters 
(10 « P « 100 d), which are often joined by smaller 
companion planets (Huang et al. 2016). Hot Jupiters 
with nearby planets may be key in understanding the 
connection of formation pathways to the observed hot 
Jupiter population. 

The Transiting Exoplanet Survey Satellite (TESS, 
Ricker 2015) is well suited to the discovery of hot 
Jupiters and potential nearby companions, as its almost- 
all-sky coverage is expected to observe nearly every 
known hot Jupiter system and discover hundreds or 
thousands more (Sullivan et al. 2015; Barclay et al. 
2018). This is particularly important since the hot 
Jupiter sample is currently heterogeneous and incom- 
plete (Yee et al. 2021). In addition, TESS has the pho- 
tometric precision to identify smaller planets down to 
~0.7 Rg (e.g. Kostov et al. 2019b; Gilbert et al. 2020; 
Silverstein et al. 2022). 

Here we present the T'ESS discovery of TOI-822.02 — 
henceforth referred to as WASP-132 c — a small planet 
associated with hot Jupiter WASP-132 b first discov- 
ered by Hellier et al. (2017). The new planet WASP-132 
c is on a 1.01 d orbit interior to the 7.13 d orbit of the 
hot Jupiter WASP-132 b. This makes the WASP-132 
system the fourth such system containing a hot Jupiter 
with a nearby small planetary companion, widening the 
sample of this rare subclass of hot Jupiters and fur- 
ther opening the possibility for comparative planetology 
both within the hot Jupiter system and between hot and 
warm Jupiter systems. 

Section 2 details the discovery of the WASP-132 c sig- 
nal in the TESS photometric data as well as the initial 
vetting efforts. Section 3 presents our refinement of the 
stellar parameters of the host star using a series of in- 
dependent models. Section 4 describes our methodology 
for modeling the full photometric light curve to obtain 
precise planetary and orbital parameters for WASP-132 
c as well as the confirmed hot Jupiter WASP-132 b. Sec- 
tion 5 presents the validation of WASP-132 c as a planet 
based on ground-based follow-up observations and sta- 
tistical validation software. Section 6 details the dynam- 
ical simulations modeling a two-planet WASP-132 sys- 
tem to probe the long-term stability of the system. Sec- 
tion 7 discusses the implications of the discovery of such 
a system in terms of hot Jupiter formation and the larger 
hot Jupiter sample. 


2. SIGNAL SEARCH AND VETTING 


WASP-132 (TOI-822, TIC 127530399) was observed 
by TESS in Sector 11 from UT April 23 to May 20, 2019 
(23.96 d) in CCD 2 of Camera 1 and in Sector 38 from 
UT April 29 to May 26, 2021 (26.34 d) in CCD 1 of Cam- 
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era 1. Data for WASP-132 were collected at 2 minute 
cadence in Sectors 11 and 38 and at 20 second cadence 
in Sector 38. The star was prioritized for high-cadence 
measurements as part of the Cycle 1 Guest Investiga- 
tor Programs G011112, G011183, and G011132 and Cy- 
cle 3 Guest Investigator Programs G03278, G03181, and 
G03106.? 

The TESS Science Processing Operation Center 
(SPOC) pipeline (Jenkins et al. 2016) processed the 
short cadence pixel data and generated the target pixel 
files (TPFs) and light curves cleaned of instrumen- 
tal systematics. The Transiting Planet Search module 
(TPS, Jenkins 2002; Jenkins et al. 2010; Jenkins et al. 
2020) of the SPOC pipeline searched the generated 2- 
minute light curves for periodic, transit-like signals for 
each TESS sector independently and jointly as a single 
light curve. TPS recovered the previously confirmed hot 
Jupiter WASP-132 b as well as a new signal at 1.01153 
d with a signal-to-noise ratio (SNR) of 10.6 in the com- 
bined data from the two sectors. This signal’s depth cor- 
responds to a planet with a radius of 2.33Rẹ when us- 
ing the stellar radius value for this target contained in 
the TESS Input Catalog (TICv8.2, Stassun et al. 2018, 
2019) to calculate the potential planet radius. The tran- 
sits of the hot Jupiter and WASP-132 c in the TESS 
data can be seen in Figure 1. 

To provide an independent recovery of this periodic 
signal, we searched the available WASP-132 TESS light 
curves with the Transit Least Squares (TLS) search al- 
gorithm (Hippke & Heller 2019). TLS utilizes analyti- 
cal transit shapes, making it more sensitive to planet 
transits than the conventional Box Least Squares (BLS, 
Kovács et al. 2002) search method and more finely tuned 
to detect small, short-period planets such as WASP-132 
c. 

Our transit search with TLS made use of the 
systematics-corrected Presearch Data Conditioning 
Simple Aperture Photometry (PDC_SAP) TESS light 
curves generated by the TESS SPOC pipeline (Smith 
et al. 2012; Stumpe et al. 2012, 2014) at the 2 minute 
and 20 second cadence for TESS Sectors 11 and 38, re- 
spectively’. We used the lightkurve Python package 
(Lightkurve Collaboration et al. 2018) to download the 
data from the Mikulski Archive for Space Telescopes 
(MAST). The light curve exhibited small amplitude stel- 
lar variability and was detrended using lightkurve’s 


2 Details of approved TESS Guest Investigator Programs 
are available from  https://heasarc.gsfc.nasa.gov/docs/tess/ 
approved-programs.html 

3 We elected to use the shortest cadence available for each sector 
to capture the shape of the transit as accurately as possible. 


built in flatten method. A window length of 70.5 days 
was chosen as it was large enough to preserve transit sig- 
nals (of duration on the order of <3 hours) while small 
enough to remove the slight stellar variability present in 
the light curve. 

Using TLS on the clean light curve, we recovered the 
known signal of hot Jupiter WASP-132 b as well as a 
signal with a period of 1.0119 + 0.0032 d with a false 
alarm probability (FAP) < 1074, which is well below the 
threshold of what Hippke & Heller (2019) states is a sig- 
nificant detection above white noise. The period, depth, 
and mid-transit time of this recovered signal are consis- 
tent with the values reported by the SPOC pipeline and 
listed on the Exoplanet Follow-up Observing Program- 
TESS (ExoFOP-TESS, ExoFOP 2019) website*. This 
recovery with TLS served as an independent check to 
ensure that the signal was not a pipeline-specific detec- 
tion. 

We performed multiple initial checks of the signal and 
TESS light curves for astrophysical false-positive scenar- 
ios that can mimic exoplanet transits. The Data Valida- 
tion module (DV, Twicken et al. 2018; Li et al. 2019) 
of the SPOC pipeline performs a suite of diagnostic vet- 
ting tests to investigate the likelihood of many of these 
false-positive scenarios. These tests include a depth test 
of the odd and even transits, a statistical bootstrap test 
that accounts for the non-white nature of the light curve 
to estimate the probability of a false alarm from random 
noise fluctuations, a ghost diagnostic test to compare the 
detection statistic of the optimal aperture against that 
of a halo with a 1 pixel buffer around the optimal aper- 
ture, and a difference image centroid test. WASP-132 
c passed all of these diagnostic vetting tests. Addition- 
ally, all TICv8.2 objects other than the target star were 
statistically excluded as sources of the 1.01 d transit sig- 
nal since the difference centroid offset tests located the 
source of the transit signal to within 1 + 3 arcsec of the 
target position. The Threshold Crossing Event (TCE) 
was promoted to TESS Object of Interest (TOI; Guer- 
rero et al. 2021) status and designated TOI-822.02 by 
the TESS Science Office based on the clean model fit 
and diagnostic test results in the SPOC data validation 
report. 

In addition to the vetting checks performed by the 
SPOC pipeline, we used DAVE (Discovery and Vetting 
of Exoplanets, Kostov et al. 2019a) to further check 
for astrophysical false-positive scenarios. DAVE is an 
automated vetting pipeline built upon many of the 
tools developed for vetting planets in Kepler data (e.g. 


4 https: //exofop.ipac.caltech.edu/tess/ 
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Figure 1. TESS PDC_SAP light curve of the WASP-132 system with the in-transit times of hot Jupiter WASP-132 b (blue) 
and WASP-132 (red) highlighted. Both Sector 11 (top) and Sector 38 (second from top) are shown. The data is detrended with 
a Gaussian Process according to the method outlined in Section 4. These Gaussian Process noise models are shown in green for 
Sector 11 (second from bottom) and Sector 38 (bottom) overlaid on the PDC_SAP light curves to show how they capture the 


variability in the light curve. 


RoboVetter, Coughlin et al. 2014). It has been used ex- 
tensively in vetting planets in K2 (Hedges et al. 2019; de 
Leon et al. 2021) and TESS (Kostov et al. 2019b; Cross- 
field et al. 2019; Gilbert et al. 2020) data as well. DAVE 
performs two sets of vetting tests. The first are light 
curve-based vetting tests searching for odd/even tran- 
sit differences, secondary eclipses, and light curve mod- 
ulations that could introduce transit-like signals. The 
second set of tests are image-based that check the pho- 
tocenter motion on the TESS image during transit. 

Unfortunately, due to the weak signal resulting 
from the shallow transit depth and relatively dim 
(Tmag-11.11) stellar host, many of the results from 
DAVE were inconclusive but still showed no significant 
indication of an astrophysical false-positive scenario. As 
such, we determined that further, more detailed model- 
ing of the transits was warranted. 


3. REFINEMENT OF HOST STAR PARAMETERS 


To perform a comprehensive analysis of the transit 
and system light curve, it was necessary to refine the 
stellar parameters of the host star, WASP-132. The 
'TICv8.2 reports key stellar parameters determined via 


independent analysis, but there are also stellar parame- 
ters reported by the WASP Collaboration in the initial 
discovery and confirmation of the hot Jupiter WASP-132 
b (Hellier et al. 2017). We performed our own indepen- 
dent analysis given available data in order to determine 
the best values for the stellar parameters to use when 
modeling the TESS data. The results of each indepen- 
dent analysis method are contained in Table 1 for com- 
parison. We find most stellar parameter values from each 
analysis method are consistent with each of the others as 
well as with both those reported by Hellier et al. (2017) 
and the TICv8.2 (Stassun et al. 2018, 2019). The excep- 
tion is a slight difference between the bolometric flux 
Foi reported by the two SED analyses. The adopted 
stellar parameter values used in the remaining validation 
and analysis of WASP-132 c are based on the isochrone 
analysis described below and are contained in Table 2. 
We also note that we see a low-amplitude 8 d peri- 
odic variation upon visual inspection of the Sector 38 
light curve that does not match up with the 33 day stel- 
lar rotation period stated in Hellier et al. (2017). If the 
8 d variability were to represent the rotation period of 


the star, this would imply a v sin i of ~ 5 km s71, us- 
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ing the equation v sin i = (27R)/P and assuming the 
star is viewed edge-on. T'his is well outside of the con- 
fidence interval of the measured value v sin à = 0.9 + 
0.8 km s7}, suggesting that this 8 d variability is not 
the stellar rotation period, although it may represent 
one of the harmonics of the true period if it is intrinsic 
to the star. However, extracted light curves using dif- 
ferent apertures do not contain this —8 day variation, 
suggesting that this shorter-scale variability may not be 
inherent to the WASP-132 system. Regardless of the ori- 
gin of this additional variability, our conclusions remain 
the same regarding the planetary nature of the 1.01 d 


transit signal. 

Overall, WASP-132 does not show significant signs 
of activity. There is the possible 8 d variability and 
reported 33 d rotation period, activity that occurs 
on much longer timescales than the orbital period of 
WASP-132 c. We also found no evidence of flares or star 
spot crossings in either the space- or ground-based data. 
In the modeling of the system's light curves (see Section 
4), a Gaussian Process was used to capture any variabil- 
ity in conjunction with transit models for each planet. 
Thus, photometric variability was modeled out while not 
diluting the transit signals (see Figure 1), in order to 
precisely measure the planet and orbital parameters. 


3.1. Isochrone Analysis 


We performed an isochrone-based analysis for WASP- 
132 using isoclassify (Huber et al. 2017; Berger et al. 
2020), which produces fundamental stellar parameters 
from a combination of input observables. We used spec- 
troscopic Teg and metallicity from the discovery paper 
(Hellier et al. 2017), Gaia Data Release 2 (Prusti et al. 
2016; Gaia Collaboration et al. 2018; Bailer-Jones et al. 
2018, DR2) parallax and coordinates, and the Two Mi- 
cron All-Sky Survey (Skrutskie et al. 2006, 2MASS) Ks 
magnitude as inputs. We also used the allsky extinction 
map detailed in Bovy et al. (2016) to estimate the photo- 
metric extinction based on the coordinates and distance 
inferred from the parallax. The best-fit values and their 
uncertainties are compiled in Table 1, and we estimate 
extinction to be Ay = 0.097 + 0.024 mag, which is con- 
sistent with the Ay = 0.093 + 0.031 mag reported in 
the TICv8.2. 


3.2. KGS SED Analysis 


As an independent determination of the basic stellar 
parameters, K.G. Stassun (KGS) performed an analysis 
of the broadband spectral energy distribution (SED) of 
the star together with the Gaia Early Data Release 3 
(Prusti et al. 2016; Brown et al. 2021, EDR3) parallax 
(with no systematic offset applied; see, e.g., Stassun & 


Torres 2021), in order to determine an empirical mea- 
surement of the stellar radius, following the procedures 
described in Stassun & Torres (2016); Stassun et al. 
(2017); Stassun et al. (2018). We pulled the the JH Kg 
magnitudes from 2MASS, the W1-W3 magnitudes from 
the Wide-field Infrared Survey Explorer (Wright et al. 
2010; Cutri et al. 2013, WISE), and the GgpGrp mag- 
nitudes from Gaia. Together, the available photometry 
spans the full stellar SED over the wavelength range 
0.4-10 um (see Figure 2). 
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Figure 2. Spectral energy distribution of WASP-132. Red 
symbols represent the observed photometric measurements 
outlined in Section 3.2, where the horizontal bars represent 
the effective width of the passband. Blue symbols are the 
model fluxes from the best-fit NextGen atmosphere model 
(black). 


We performed a fit using NextGen stellar atmosphere 
models, with the free parameters being the effective tem- 
perature (Ter) and metallicity ([Fe/H]), as well as the 
extinction Ay, which we fixed at zero due to the prox- 
imity of the system to Earth. The resulting fit (Fig- 
ure 2) has a best-fit Tyg = 4750 + 75 K, [Fe/H] = 
0.0 + 0.5, with a reduced x? of 0.8. Integrating the 
(unreddened) model SED gives the bolometric flux at 
Earth, Fyo1 = 5.442+0.063 x 107 !? erg s^! cm ?. Taking 
the Fyo; and Teg together with the Gaia parallax gives 
the stellar radius, R, = 0.752 + 0.024 Ro. In addition, 
we can estimate the stellar mass from the empirical rela- 
tions of Torres et al. (2010), giving M, = 0.80-£0.05 Mo. 
Finally, the reported stellar rotation period of ~33 d 
implies an age of 7, = 3.2 + 0.5 Gyr via the empiri- 
cal gyrochronology relations of Mamajek & Hillenbrand 
(2008). 


3.3. MLS SED Analysis 
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Table 1. Stellar parameters obtained using each of the methods outlined in Section 3. Values from Hellier et al. (2017), which discovered 
WASP-132 b, and from the TICv8.2 (Stassun et al. 2018, 2019) are included for comparison. The final adopted parameters for this 
analysis are contained in Table 2. All uncertainties reported are the lo value. 


Parameter 'TICv8.2 Hellier et al. (2017) Isochrone Analysis KGS SED Analysis MLS SED Analysis 
Test (K) 4742 + 129 4750 + 100 4714187 4750 + 75 4753 + 80 
[Fe/H] — 0.22 + 0.13 0.18 + 0.12 0.0 + 0.5 — 
M. (Mo) 0.760 + 0.089 0.80 + 0.04 0.782 + 0.034 0.80 + 0.05 — 
Rs (Ro) 0.790 + 0.057 0.74 + 0.02 0.753* 0:028 0.752 + 0.024 0.767 + 0.026 
L, (Lo) 0.284 + 0.012 — 0.253* 9032 — 0.271 0.007 
log(g) 4.524 0.094 4.6 + 0.1 4.5767 O36 — — 
px (g cm?) 2.17 + 0.59 2821020 1811545 — — 
Age (Gyr) = >0.5 7.055 * 2-8 3.2 + 0.5 = 
Distance (pc) 122.91 + 0.57 120 + 20 126 +5 — — 
Foo (erg s! cem ? x10 19) — — — 5.442 + 0.063 5.69 + 0.14 


As an additional check on our stellar effective temper- 
ature, luminosity, and radius results, M. L. Silverstein 
(MLS) led a second SED analysis following methodol- 
ogy similar to Dieterich et al. (2014) (Silverstein et al., 
in preparation). Spanning optical to mid-infrared wave- 
lengths, we extract JH Ks magnitudes from 2MASS and 
W1-W3 magnitudes from WISE as in the previous sub- 
section. We differ in our adoption of the Gaia DR2 par- 
allax and of V RI photometry converted from Gaia DR2 
GaGgpGgp. We compare nine different color combina- 
tions to the B T-Settl 2011 photospheric models (Allard 
et al. 2012) to derive a best-fitting Tag = 4753 + 80 K, 
assuming [Fe/H] = 0. Next we iteratively scale the 
resulting best-match model using a polynomial func- 
tion until model and observed photometry match to 
within the error bars. We then integrate the final spec- 
trum and apply a bolometric correction to determine 
Fyo1 = 5.691 + 0.137 x 1071? erg s^! em~?, and we scale 
by the parallax to derive L, = 0.271 + 0.007 Lo. A ra- 
dius of R, = 0.767 + 0.026 Re is then calculated using 
the Stefan-Boltzmann Law. These results are listed in 
Table 1 and match those from the other independent 
methods described in this paper. 


4. MODELING THE PHYSICAL PROPERTIES OF 
WASP-132 C 


While TLS is useful at detecting signals, the period 
grid that it searches is not very fine by default. Com- 
bined with the refined stellar parameters discussed in 
Section 3, it is possible to model the light curve in a 
more detailed fashion than the initial transit search to 
find the maximum likelihood values for the planet prop- 
erties in the system. To perform this detailed modeling, 
we used the software exoplanet (Foreman-Mackey et al. 


2019). exoplanet is a toolkit for probabilistic modeling 
of transit and radial velocity observations of exoplanets 
using PyMC3. This is a powerful and flexible program 
that can be used to build high-performance transit mod- 
els and then sample them through Markov Chain Monte 
Carlo (MCMC) simulations to provide precise transit 
and orbital parameters. 

We utilized the same PDC_SAP light curves used in 
the transit search with TLS with one difference. For the 
modeling with exoplanet, we did not apply any initial 
detrending that could possibly alter the transit signals 
but instead included a Gaussian Process (GP) in the 
model. Our model had three elements: two planet com- 
ponents with Keplerian orbits and limb-darkened tran- 
sits (one for each potential planet in the system) and 
a GP component that modeled residual stellar variabil- 
ity. The planet models were computed using STARRY 
(Luger et al. 2019) and the GP was computed us- 
ing celerite (Foreman-Mackey et al. 2017; Foreman- 
Mackey 2018). The GP component models the resid- 
ual stellar variability in the light curve and describes a 
stochastically-driven, damped harmonic oscillator with 
two hyper-parameters, In(So) and In(wo), which repre- 
sent the undamped angular frequency of a simple har- 
monic oscillator and the power at w — 0, respectively. 
We fixed the quality factor Q of the simple harmonic 
oscillator to 1 / V2 and put wide Gaussian priors on 
In(Sg) and In(wo), setting their means to the natural 
log of the standard deviation of the flux and natural log 
of one tenth of a cycle, respectively, with both of their 
standard deviations set to 10. This form of GP has the 
advantage of being able to model a wide range of low 
frequency astrophysical and instrumental signals with- 
out requiring a physical model for the observed variabil- 
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'Table 2. Adopted stellar parameters for WASP-132. 


Parameter Value Source 
Identifying Information 
Name WASP-132 
TIC ID 127530399 'TICv8.2 
TOIID TOI-822 Guerrero et al. (2021) 
Alt. Name UCACA 220-083803 


a R.A. (hh:mm:ss) 
ó Dec. (dd:mm:ss) 


jia (mas yr?) 
is (mas yr) 
Distance (pc) 


Spectral Type 
Tere (K) 


Rotation period (d) 


v sin i (km s!) 
Age (Gyr) 


B (mag) 
V (mag) 
Ga (mag) 
Gerp (mag) 


Astrometric Properties 
14:30:26.21 (J2015.5) 
-46:09:34.29 (J2015.5) 


12.255 + 0.020 
-73.169 + 0.022 
123.57 + 0.29 


Stellar Properties 
K4 
4714+87 
0.18 + 0.12 
0.782 + 0.034 
0.753* 9:626 
0.253* 9022 
4.57679 036 
1.81*21$ 
7.33 
0.9 + 0.8 
3.2 + 0.5 
Photometric Properties 
13.142 + 0.011 
11.938 + 0.046 
11.7467 + 0.0002 
12.3000 + 0.0007 
11.0487 + 0.0004 


11.111 + 0.006 
10.257 + 0.026 
9.745 + 0.023 
9.674 + 0.024 
9.557 + 0.022 
9.638 + 0.020 
9.575 + 0.040 
8.281! 


Gaia EDR3 
Gaia EDR3 
Gaia EDR3 
Gaia EDR3 
Gaia EDR3 


Hellier et al. (2017) 
This Work 
This Work 
This Work 
This Work 
This Work 
This Work 
This Work 
Hellier et al. (2017) 
Hellier et al. (2017) 
This Work 


APASS DR9 
APASS DR9 
Gaia EDR3 
Gaia EDR3 
Gaia EDR3 
TICv8.2 
2MASS 
2MASS 
2MASS 
AIIWISE 
AIIWISE 
AIIWISE 
AIIWISE 


Gaia EDR3 - Prusti et al. (2016); Brown et al. (2021), TICV8.2 - 
Stassun et al. (2019), APASS DR9 - Henden et al. (2016), 2MASS - 
Skrutskie et al. (2006), AIIWISE - Cutri et al. (2013). !Only a limit 
is reported for W4 since the signal-to-noise ratio was too low for a 
confident detection. 


ity. We also included a white noise term in the model 
which is parameterized by the natural log of the stan- 
dard deviation of the flux with a prior identical to that 
of In(So). The GP parameters of the two TESS sectors 
were modeled separately since the sectors may have dif- 
ferent noise parameters, especially since the data were 
taken at different cadences. 

'The planet model was parameterized with a two term 
limb-darkening component and the stellar radius and 
mass. The individual planets were parameterized in 
terms of the natural log of orbital period, mid-transit 
time, transit depth, impact parameter, eccentricity, and 
periastron angle at time of transit. For our priors on 
the stellar parameter components of the model, we used 
the mean and standard deviation values of our anal- 
ysis discussed in Section 3 and displayed in Table 2. 
We followed Kipping (2013a) for the parameterization 
of the limb-darkening. We used the SPOC values listed 
on ExoFOP-TESS for the means and standard devia- 
tions of Gaussian priors on the natural log of the orbital 
period, mid-transit time, and transit depth for the two 
planets. We imposed a uniform prior on the impact pa- 
rameter bounded between 0 and 1. For the eccentricity 
prior, we used a Beta prior with a = 0.867 and 8 = 3.03 
as suggested by Kipping (2013b). The eccentricity was 
bounded between zero and one and sampled as ecos(w). 
The periastron angle at transit was sampled in vector 
space to avoid the sampler seeing discontinuities. We 
sampled the posterior distribution of the model parame- 
ters using the No U-turn Sampler (NUTS, Hoffman et al. 
2014), which is a form of Hamiltonian Monte Carlo, as 
implemented by PyMC3 (Salvatier et al. 2016a). We ran 
3 simultaneous chains, each with 2000 tuning steps and 
2500 draws in the final sample. 

Initially, since individual TESS sectors often have dif- 
ferent noise properties, we modeled both Sectors 11 and 
38 independently from one another using the model de- 
scribed above. However, the resulting posterior distri- 
butions were equivalent within lo errors, so we decided 
to combine both sectors of TESS data into a single light 
curve and use the same model. Since the two light curves 
are separated by >1 year, modeling both together as a 
single light curve increases the time baseline with which 
to model the orbit of the system, allowing for a better 
constrained orbital period than any individual sector or 
two sectors back-to-back in time. We binned the Sec- 
tor 38 TESS data from 20 second cadence to 2 minute 
cadence to match with the TESS Sector 11 data in 
order to create a uniform data set with an increased 
photometric precision in Sector 38. No correction for 
contamination from nearby sources of constant bright- 
ness was included since follow up observations by SOAR 
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Table 3. Planet and orbital parameters for 
WASP-132 b and c calculated by modeling the 
TESS photometric data with exoplanet. Er- 
rors are reported at the lo level. Noise param- 
eters are also included. 


Parameter Value 


Model Parameters 


Star 

Limb darkening ui 
Limb darkening u2 
Radius [Ro ] 
Mass [Mo ] 
In(pap) 

In(o GP,S11) 
In(cGp,ssa) 


WASP-132 c 

To (BJD - 2457000) 
In(Period) [days] 
Impact parameter 
In( Transit Depth) 
eccentricity 


w [radians] 


WASP-132 b 

To (BJD - 2457000) 
In(Period) [days] 
Impact parameter 
In(Transit Depth) 
eccentricity 


w [radians] 


0.43 + 0.11 
0.17 + 0.24 
0.754 + 0.024 
0.781 + 0.033 
1.03 + 0.15 
-7.96 + 0.16 
-7.11 + 0.16 


1597.5762 + 0.0024 
0.011 + 4.69e-6 
0.28*0-1$ 
-7.437 + 0.068 
0.13*0:28 
—0.93*1-82 


2337.6080 + 0.0002 


1.96 + 5.49e-7 
0.16 + 0.11 
-4.026 + 0.01 


0.07. 0-8 
—0.08*2-55 


Derived Parameters 


WASP-132 c 
Period [days] 
Ry/ R. 

Radius [Re | 
Radius [Ry] 
a/R. 

a [AU] 
Inclination [deg] 


Duration [hours] 


WASP-132 b 
Period [days] 
Ry/ R. 

Radius [Ra | 
Radius [Ry] 
a/R. 

a [AU] 
Inclination [deg] 


Duration [hours] 


1.011534 + 0.000005 
0.023 + 0.001 
1.85 + 0.10 
0.165 + 0.009 
5.17 + 0.18 
0.0182 + 0.0003 


86.64*112 
1.470: 


1.133514 + 0.000004 
0.122 + 0.006 
10.05 + 0.34 
0.897 + 0.030 
19.03 + 0.66 
0.067 + 0.001 


89.51*0-18 
3.187035 


Normalized Flux 


-0.15 —0.10 —0.05 0.00 0.05 0.10 0.15 
Time Since Transit [days] 


Figure 3. TESS data (light gray) for WASP-132 c (top) and 
WASP-132 b (bottom) phase folded to the best-fit period and 
mid-transit time with the exoplanet model overlaid. The 
process of fitting the transit model is described in Section 
4. The blue points are the phase-folded photometric data 
binned for clarity. 


and LCOGT cleared the nearby field of contaminating 
sources of this nature (see Section 5.1). 

The median values and 1c errors for the best-fit tran- 
sit model are contained in Table 3 and the folded light 
curve with the best-fit and lo bounds of the transit 
model are shown in Figure 3. 


5. VALIDATION OF WASP-132 C 


While the TESS pipeline and DAVE perform vetting 
analysis against possible false alarm and false positive 
scenarios, they alone are not rigorous enough to validate 
the planetary nature of WASP-132 c. We therefore inves- 
tigated this signal using both observational constraints 
(Section 5.1) as well as publicly-available statistical soft- 
ware packages vespa and TRICERATOPS (Section 5.2). 


5.1. Follow-up Observations 


In order to better constrain the false positive proba- 
bility of a planet candidate, follow-up observations can 
rule out sections of the parameter space of different as- 
trophysical false positive scenarios. Since WASP-132 is 
already known to host a confirmed hot Jupiter, both 
speckle imaging and radial velocity follow-up observa- 
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TIC127530399 
1 SOAR Speckle ACF 


Aó [arcsec] 


Amagnitude (/-band) 


0.0 0.5 1.0 1.5 2.0 2.5 3.0 
Aarcsec 


Figure 4. The 5-c detection sensitivity of the SOAR speckle 
imaging of WASP-132, with inset two-dimensional auto- 
correlation function reconstructed image of the field. The 
data indicate that there are no close-in companions within 3 
arcsecond of WASP-132. 


tions were readily available for our analysis and could 
be included in our final validation of WASP-132 c. 


5.1.1. SOAR Speckle Imaging 


High-angular resolution imaging is needed to search 
for nearby sources that can contaminate the TESS pho- 
tometry (resulting in an underestimated planetary ra- 
dius) or be the source of astrophysical false positives, 
such as background eclipsing binaries. We searched for 
stellar companions to WASP-132 with speckle imaging 
using the HRCam instrument on the 4.1-m Southern 
Astrophysical Research (SOAR) telescope (Tokovinin 
2018) on 10 February 2020 UT, observing in Cousins 
I-band, a similar visible bandpass as TESS. This obser- 
vation was sensitive to a 5.0-magnitude fainter star at 
an angular distance of 1 arcsec from the target. More 
details of the observation are available in Ziegler et al. 
(2020). The 5c detection sensitivity and speckle auto- 
correlation functions from the observations are shown 
in Figure 4. No nearby stars were detected within 3” of 
WASP-132 in the SOAR observations. 


5.1.2. CORALIE Radial Velocity Observations 


In addition to SOAR speckle imaging, prior to the dis- 
covery of WASP-132 c there already existed radial ve- 
locity measurements of the host star WASP-132 which 
were used to confirm the planetary nature and mea- 
sure the mass of the hot Jupiter WASP-132 b. In total, 
there were 36 radial velocity measurements taken across 
a time span of almost 2 years. All of these measurements 


were obtained with the CORALIE spectrograph which 
is an echelle spectrograph mounted on the 1.2-m Euler 
telescope in La Silla, Chile. These data are published in 
Hellier et al. (2017), which provides further information 
regarding the method used to reduce the radial velocity 
data. 

According to the forecaster Python package (Chen 
& Kipping 2016), the projected mass of WASP-132 c 
should be 4.45 Ma. forecaster assumed a gaseous en- 
velope when estimating the mass of WASP-132 c, how- 
ever the radius of the planet may place it in the super- 
Earth regime of planets. This would imply that the den- 
sity of the planet is higher than a gaseous planet, mean- 
ing that the mass follows the relation My ~ a. This 
would result in a mass of 9.74 Mẹ. We made no dis- 
tinction between these two scenarios in our analysis of 
CORALIE data since we imposed a very wide prior on 
the mass of WASP-132 c. 

We first performed a joint fit of the time series TESS 
photometry with the CORALIE radial velocity measure- 
ments using exoplanet. This way, we were able to take 
advantage of the unique strengths of each data type in 
constraining the orbital and planet parameters as well as 
accounting for the noise and variability in the data. We 
used a model similar to that described in Section 4 but 
with the addition of planet mass as well as a quadratic 
trend and a jitter term for the radial velocities. The pri- 
ors on the trend and jitter model components were Nor- 
mal distributions with the quadratic trend distribution 
centered on 0 and the jitter term distribution centered 
on the standard deviation of the radial velocity data. We 
imposed a wide log-uniform prior on the mass from 0 to 
8. This is more than wide enough to encompass the pro- 
jected masses of both potential compositions for WASP- 
132 c as well as the 30 upper limit on the reported mass 
of WASP-132 b. CORALIE also underwent an upgrade 
partway through the dataset and this is accounted for by 
an offset between the pre- and post-upgrade portions of 
the data. We also fit for the trend in the data described 
in Hellier et al. (2017). 

The posterior probability distribution for the mass of 
WASP-132 c was highest around zero with a wide spread 
of values, indicating a non-detection in the radial ve- 
locity. Therefore, we report a 3e upper limit from this 
distribution of 37.35 Ma, corresponding to a radial ve- 
locity semi-amplitude of 28.23 ms~!. All other planetary 
and orbital parameters modeled with the joint TESS 4- 
CORALIE data (e.g. period, mid-transit time, radius, 
etc.) are consistent with the values obtained from mod- 
eling the TESS data alone, described in Section 4. The 
best fit models are plotted with the phase-folded radial 
velocity data for WASP-132 b and c in Figure 5. We note 
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Figure 5. The CORALIE radial velocity measurements 
phase-folded to the best-fit periods of WASP-132 c (top) and 
WASP-132 b (bottom). The solid orange line represents the 
median radial velocity model and the shaded regions repre- 
sent the lo uncertainties in the model. 


that Hellier et al. (2017) report that the CORALIE ra- 
dial velocities show excess scatter, which may be due 
to magnetic activity. Therefore, the upper limit on the 
mass reported here may be inflated due to the high scat- 
ter. 

As a check on the validity of the procedure, we also 
compared the measured planet mass for the confirmed 
hot Jupiter in the system WASP-132 b against its re- 
ported mass in Hellier et al. (2017). We measure the 
mass to be 121.89 + 21.85 Mẹ which is in agreement 
with the 130.31 + 9.53 Mẹ measured by Hellier et al. 
(2017). The uncertainty on our calculated mass of the 
hot Jupiter is higher than that reported by Hellier et al. 
(2017) likely because Hellier et al. (2017) also joint fit 
the CORALIE data with photometric data from both 
Wide Angle Search for Planets (WASP) and TRAnsit- 
ing Planets and PlanetesImals Small Telescope (TRAP- 
PIST) observations that our analysis does not include 
as well as other ancillary RV observation data products. 

As an independent check of the values obtained from 
exoplanet, we used the online Data and Analysis Cen- 


ter for Exoplanets? (DACE) platform to model the ra- 
dial velocity data and investigate the significance of any 
signals around 1.01 d (the period of WASP-132 c). The 
Keplerian model initial conditions for the hot Jupiter 
were input based on the values on ExoFOP-TESS for 
WASP-132 b and a decrease of 60 ms~! over the course 
of the observations was added as noted in Hellier et al. 
(2017). When viewing the periodigram of the radial ve- 
locity data after adding in these components, there ap- 
pear to be spikes in signal around 1 d, but none of which 
has an FAP « 1096. This is not surprising, since many 
signals are aliased to 1 d and this is à common phe- 
nomenon in analyzing radial velocity data. The addition 
of another Keplerian orbit at the expected orbital period 
of WASP-132 c with the predicted semi-amplitude based 
on the mass predicted by forecaster results in a higher 
reduced x?, indicating a worse fit. 

According to Equation 14 in Lovis & Fischer (2011), 
the expected semi-amplitude of the radial velocity sig- 
nal of WASP-132 c with 4.45 Ma as predicted by 
forecaster is 3.36 ms-!. If an Earth-like composition 
is assumed, the mass of 9.74 Mẹ would result in a semi- 
amplitude of 7.44 ms~'. CORALIE has been reported 
to have an individual measurement precision ranging be- 
tween 3.5 and 6 ms ^! (Rickman et al. 2019), putting 
the projected mass of WASP-132 c slightly below and 
slightly above the lower limit of CORALIE's detectable 
parameter space for gaseous and Earth-like composi- 
tions, respectively. Combined with the fact that the or- 
bital period falls very close to the highly-aliased value of 
1 d, it is logical that there was no significant detection 
of WASP-132 c in the CORALIE radial velocity data. 
Further radial velocity observations are necessary with a 
more precise instrument such as the High Accuracy Ra- 
dial velocity Planet Searcher (HARPS) or the Carnegie 
Planet Finder Spectrograph (PFS) in order to obtain a 
mass measurement for the planet WASP-132 c. 


5.1.8. LCOGT 1m 


We observed T'OI-822 c from the Las Cumbres Obser- 
vatory Global Telescope (LCOGT; Brown et al. 2013) 
1.0m network node at the South Africa Astronomical 
Observatory on UTC 2022 March 5, 2022 March 8, 
and 2022 March 9 in Sloan 7’ band. We used the TESS 
Transit Finder, which is a customized version of the 
Tapir software package (Jensen 2013), to schedule our 
transit observations. The 1m telescopes are equipped 
with 4096 x 4096 SINISTRO cameras having an image 
scale of 07389 per pixel, resulting in a 26’ x 26' field 
of view. The images were calibrated by the standard 


5 https://dace.unige.ch 
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LCOGT BANZAI pipeline (McCully et al. 2018). Photo- 
metric data were extracted using AstroImageJ (Collins 
et al. 2017). The target star apertures exclude virtu- 
ally all of the flux from the nearest Gaia EDR3 and 
TESS Input Catalog known neighbor (TIC 1051778874), 
which is 10"0 northwest of TOI-822. The light curve data 
are available on the ExoFOP-TESS website? (ExoFOP 
2019). 

'The UTC 2022 March 5 observation used focused ob- 
servations that were intended to saturate the target star 
for purposes of checking fainter stars within 2/5 of TOI- 
822 c for a potential nearby eclipsing binary (NEB) that 
could be the cause of the periodic detection in the T'ESS 
data. We rule out NEB signals, relative to the depth re- 
quired in each neighboring star given its TIC version 8 
delta TESS magnitude, by a factor of more than three 
times the RMS of the light curve scatter for all neigh- 
boring stars out to 50" from TOI-822. This is consistent 
with the SPOC pipeline’s (lack of) centroid offset find- 
ing that the source of the TESS-detected signal is within 
~ 10” of the target star location using data from TESS 
sectors 11 and 38. We also do not see any obvious NEB 
signals in stars out to 2/5, although some light curves 
suffer from blending from brighter nearby stars. 

Although we expected the target star to be saturated 
in the UTC 2022 March 5 observation, it ultimately 
was strongly exposed, but not saturated. While the de- 
trended light curve is consistent with a 40 min (2.70) 
late ~ 600 ppm deep event on-target relative to the 
SPOC sectors 11 and 38 ephemeris, the detection is con- 
sidered inconclusive due to limited post-transit baseline 
coverage and apparent systematics at the level of ~ 500 
ppm in the undetrended light curve. 

The UTC 2022 March 8 and UTC 2022 March 9 ob- 
servations were intentionally defocused to attempt to 
confirm the tentative 2.70 late ~ 600 ppm deep event in 
the UTC 2022 March 5 data. However, these data also 
suffered from ~ 500 ppm systematics in the undetrended 
light curves. With our best detrending efforts, the UTC 
2022 March 8 would marginally suggest a roughly 100 
min late ~ 600 — 700 ppm deep ingress at the end of the 
light curve. On the other hand, the UTC 2022 March 9 
detrended light curve would suggest a roughly 20 min 
early ~ 700 ppm deep event. 

Given the level of systematics in the undetrended data 
and the inconsistent timing of the three tentative tran- 
sit detections in the detrended data, we do not further 
consider the on-target results in the following analyses. 
Although we favor the interpretation that the extracted 


6 https:/ /exofop.ipac.caltech.edu/tess/target.php?id—127530399 


events are systematics driven, we cannot rule out the 
interpretation that some or all of the tentative detec- 
tions are astrophysical and that the timing offsets are 
indicative of large T TVs in the system. 


5.2. Validation Using Software Tools 


While observational constraints can rule out portions 
of the parameter space where astrophysical false posi- 
tives could exist, these observational limits are incom- 
plete and do not rule out the entirety of the false positive 
parameter space. We are able to statistically analyze 
the remaining likelihood of false positive signals using 
publicly-available software. Using the available follow- 
up observations as constraints, we ran vespa (Morton 
2012, 2015) and TRICERATOPS (Giacalone & Dressing 
2020; Giacalone et al. 2021) to further establish the plan- 
etary nature of this signal. 

vespa compares the transit signal to a number of as- 
trophysical false-positives including an unblended eclips- 
ing binary (EB), a blended background EB (BEB), 
a hierarchical companion EB, and EB scenarios with 
a double-period. We ran vespa using the TESS light 
curves detrended with the noise model described in Sec- 
tion 4 and phase-folded using the median values of the 
posteriors for the period and mid-transit time. The mod- 
eled transits of hot Jupiter WASP-132 b were also sub- 
tracted from the light curve. We included the maximum 
possible secondary depth of phase-folded features calcu- 
lated by DAVE (Section 2) and the SOAR I-band contrast 
curve (Section 5.1.1) as observational constraints when 
calculating the false positive probability of the WASP- 
132c signal. By default, vespa simulates the background 
starfield within 1 square degree of the target, but we 
dictated that the maximum aperture radius interior to 
which the signal must be produced be 42", which is the 
radius of two T'ESS pixels and the maximum size of the 
aperture that the SPOC pipeline used to extract the 
light curves. 

Using these inputs, we calculated the false positive 
probability (FPP) of the WASP-132 c signal to be 
0.00193. The only false positive scenario with any re- 
maining probability was the case of a blended back- 
ground EB, however the probability for that scenario 
was «0.01 and is highly disfavored over the planet sce- 
nario. Since the overall FPP «0.01, this signal can be 
considered statistically validated by vespa. 

As an independent check, we ran the WASP-132 
c signal through the statistical validation software 
TRICERATOPS. This software is similar to vespa in that 
it checks the signal against a set of scenarios that could 
produce transit-like signals. TRICERATOPS was specifi- 
cally designed for TESS observations and accounts for 
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known nearby stars contained within the light curve ex- 
traction aperture as well as stars within 2.5' of the tar- 
get. This tool calculates both the FPP of the signal and 
the probability that the planet candidate is a false pos- 
itive originating from a known nearby star, labeled the 
nearby FPP (NFPP). 

Since TRICERATOPS is sensitive to the extraction aper- 
ture, we ran TRICERATOPS using light curves extracted 
with two separate apertures for comparison and qual- 
ity check. The first run used the apertures that the 
TESS SPOC pipeline used to extract the PDC. SAP 
light curves, which were 12 and 10 pixels centered on 
the target position for Sectors 11 and 38, respectively. 
'The second run used custom light curves extracted using 
reduced apertures that were smaller than those used by 
the TESS SPOC pipeline with 4 and 6 pixels centered on 
the target position for Sectors 11 and 38, respectively. 

For each TRICERATOPS run, we input the apertures 
used to extract the light curves, the I-band contrast 
curve from SOAR, and the median values for period 
and mid-transit time from the modeling performed in 
Section 4. The light curves used were phase-folded at 
these period and mid-transit time values. For each set 
of apertures, we ran TRICERATOPS twenty times and took 
the average FPP values. Using the apertures from the 
TESS SPOC pipeline, we obtained an FPP = 0.0126 + 
0.0003 for WASP-132 c. Using the reduced apertures, we 
obtained an FPP = 0.0107 + 0.0004. Since the nearby 
field has been cleared of nearby and background eclips- 
ing binary systems by the LCOGT observations (see Sec- 
tion 5.1.3), SOAR, observations (see Section 5.1.1), and 
the the TESS SPOC centroid offset tests (which put the 
signal within 1 + 3 arcsec of the target), we assume an 
NFPP value of 0. 

According to Giacalone et al. (2021), for a signal to 
be statistically validated as a planet by TRICERATOPS, it 
must have FPP < 0.015. Our signal of 1.01 d meets the 
TRICERATOPS FPP criterion for both sets of extraction 
apertures. For consistency, we also reran vespa on the 
reduced aperture light curve obtained an FPP = 9.02 x 
1075, still well below the vespa statistical validation 
threshold. This difference in FPP between vespa and 
TRICERATOPS may stem from the fact that TRICERATOPS 
uses the actual TESS aperture and background star pop- 
ulation while vespa simulates this with a extraction ra- 
dius and TRILEGAL simulations. Therefore, given the 
vespa validation, sufficiently low TRICERATOPS FPP val- 
ues, and strong constraints on nearby background stars 
in the field of view, we consider this signal to be statis- 
tically validated as a planet. 

Furthermore, neither of these statistical validation 
packages account for the fact that this signal is a part 


of a multi-planet system. Lissauer et al. (2012, 2014) 
demonstrated that false positives are less likely in multi- 
planet systems and that the FPPs calculated without 
accounting for this fact should be treated as upper lim- 
its. An analysis of TESS multiplanet systems indicates 
that this “multiplicity boost” may reduce these FPPs 
by a factor of ~ 20x (Guerrero et al. 2021), although 
the scarcity of inner companions to hot Jupiters suggests 
that the multiplicity boost should be smaller than this 
factor. 


6. DYNAMICAL STABILITY 


To probe the dynamical stability of the system, we 
simulated the system using REBOUND (Rein & Liu 2012), 
a flexible N-body integrator written in both Python and 
C. We ran 30 iterations of the system, each time per- 
turbing the mass, eccentricity, inclination, and periods 
of each of the planets within their 30 uncertainties to 
test stability at the extremes of the possible parameter 
space for both planets. Each iteration was integrated for 
simulation timescales of 100 Myr and a timestep of 0.2 
d (~20% of WASP-132 c’s best-fit orbital period) using 
the MERCURIUS integrator (Rein et al. 2019). 

The MERCURIUS integrator is a hybrid integrator 
that uses the symplectic Wisdom-Holman integrator 
WHFast when particles are far apart from each other 
and switches to the higher order integrator IAS15 dur- 
ing close encounters, which integrates with a smaller, 
adaptive timestep. We specified the minimum timestep 
for the IAS15 integrator as 0.02 d in order to speed up 
computation time. 

We found that the planetary system is stable for all 
portions of the parameter space simulated on timescales 
of 100 Myr. None of the simulations exhibited a drastic 
change in semi-major axis outside of the normal oscilla- 
tions due to gravitational interactions between the three 
bodies in the system (the two planets and the host star). 
Furthermore, no collisions were registered between any 
of the bodies over the course of the simulations. 

We also simulated the system using the 3c upper lim- 
its on mass, eccentricity, and inclination for both of the 
planets in order to maximize the likelihood of a close 
encounter. We integrated this system using the MER- 
CURIUS integrator for 1 Myr at a timestep of 1.2 hours 
(~5% of WASP-132 c’s best-fit orbital period). The po- 
sitions of both planets at 1000 different evenly-spaced 
timesteps can be seen in Figure 6. We found that the 
system was stable at the extreme end of the parameter 
space of both planets as there were no close encounters 
between the planets and no significant change in the 
semi-major axis the two planets over the course of the 
simulation. This is illustrated by the gap between the 
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Figure 6. Position of WASP-132 b and c at 1000 evenly spaced points in time during the 1 Myr simulation with mass, 
eccentricity, and inclination initial values set to the 3o upper limits for both planets. WASP-132 b is denoted by the outer, 
green points and WASP-132 c is denoted by the inner, orange points. Left: A side-on view of the system as it would be seen 
from Earth. Right: A top-down view to illustrate the positions of both planets as they orbit over the course of 1 Myr. Note the 


gap between the maximum distance of WASP-132 c and the minimum distance of WASP-132 b from the host star. 


innermost positions of WASP-132 b and outermost po- 
sitions of WASP-132 c in Figure 6. Given the stability 
from both the simulations with random draws of param- 
eters as well as the simulation with the planet parameter 
values most likely to cause a close encounter, the addi- 
tion of WASP-132 c into the WASP-132 system does not 
appear to affect the long-term stability of the system. 


7. DISCUSSION 


The discovery and validation of WASP-132 c places 
the WASP-132 system among only a handful of sys- 
tems with nearby companions to a hot Jupiter. Prior to 
WASP-132, only WASP-47 (Becker et al. 2015), Kepler- 
730 (Canas et al. 2019), and TOI-1130 (Huang et al. 
2020) were known to harbor a hot Jupiter with nearby 
companions. Figure 7 illustrates all of the currently 
known systems containing hot Jupiters with nearby 
companions. These four multi-planet systems are still 
consistent with occurrence rate estimates provided by 
both Huang et al. (2016), Zhu & Dong (2021), and Hord 
et al. (2021) which calculate 1.1*1?396, ~2% (<9.7%, 
95% upper limit), and 7.3* 2296 of hot Jupiters to have 
nearby companions, respectively. 


7.1. Potential Formation Pathways 
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Figure 7. À schematic depiction of the four known systems 
with a hot Jupiter and at least one nearby companion planet 
(WASP-47, Kepler-730, TOI-1130, WASP-132). Orbital pe- 
riods of the hot Jupiters increase from top to bottom. Sizes 
of circles for the planets are to scale with one another but not 
with distance from host star and host star radius. Likewise, 
circles for the host stars are to scale with one another but not 
with planets or orbital distance. Hot Jupiters are denoted by 
red circles, companion planets by gray circles, and the stellar 
hosts are colored based on their surface temperature. 


The existence of nearby companions in four hot 
Jupiter systems suggests a dynamically cooler formation 
mechanism than the high-eccentricity migration that is 
invoked to explain a significant fraction of hot Jupiters. 


Terr [K] 
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Although TESS appears to be improving our under- 
standing of the architecture of hot Jupiter systems — 
now with a second hot Jupiter companion discovery pre- 
sented here — hot Jupiters with companions are still a 
rarity, comprising only a small percentage of the nearly 
500 hot Jupiters currently confirmed. Few hot Jupiter 
systems have yet to been searched for companions to 
this level of sensitivity, though. This scarcity of hot 
Jupiters with nearby companions suggests that high- 
eccentricity formation scenarios may dominate the ob- 
served hot Jupiter population. However, there is increas- 
ing evidence that there must be a subpopulation of hot 
Jupiters that do not form via high-eccentricity migra- 
tion. This is evidenced not only by the existence of hot 
Jupiters with nearby companions as presented here, but 
also statistical simulations and analytical calculations 
show that high-eccentricity migration alone cannot re- 
produce the observed hot Jupiter population and require 
a mechanism such as disk migration (Dawson et al. 2014; 
Anderson et al. 2016; Munoz et al. 2016). Although disk 
migration could explain these handful of unique systems, 
it is also possible that a super-Earth could have managed 
to exceed the threshold mass for runaway gas accretion, 
forming the hot Jupiter and leaving intact any planets 
orbiting interior (Lee et al. 2014; Batygin et al. 2016; 
Huang et al. 2020). 

'The subpopulation of hot Jupiters with companions 
that apparently did not undergo high-eccentricity mi- 
gration may have more in common with warm Jupiters 
(10 « P « 100 d) than their fellow hot Jupiters in terms 
of formation. In contrast to hot Jupiters, ~50% of warm 
Jupiters have nearby companion planets (Huang et al. 
2016). Hot Jupiters with companions form a seemingly 
continuous period distribution with their warm Jupiter 
counterparts (Huang et al. 2020) and may represent 
the tail end of a wider population spanning across the 
physically-unmotivated 10 d dividing line. 

Taking, for instance, a definition of hot Jupiter based 
on equilibrium temperature where Teq 71000 K defines 
a hot Jupiter (e.g. Miller & Fortney 2011; Thorngren 
& Fortney 2018) leaves only WASP-47 b and Kepler- 
730 b classified as hot Jupiters with Teq of 1259 and 
1219 K, respectively, with TOI-1130 b and WASP-132 
b classified in a slightly lower “warm” range with Teq 
of 637 and 763 K, respectively. T'his definition of hot 
Jupiters based on equilibrium temperature is equally ar- 
bitrary, however, as it still fails to provide a clear distinc- 
tion between hot Jupiter subpopulations that formed via 
different formation mechanisms. Systems from a poten- 
tial ^warm Jupiter tail", such as those discussed here, 
would inadvertently be included in this definition of hot 
Jupiters. 


7.2. Investigating System Architectures 


Comparing these unique systems to the overall exo- 
planet population (see Figure 8), it becomes apparent 
that three of the hot Jupiters with nearby companions 
have orbital periods that are longer than the typical hot 
Jupiter. The hot Jupiters with companions reside in a 
more sparsely-populated region of period-radius space 
farther from their host stars than most hot Jupiters. 
While this is notable, the sample of known hot Jupiters 
with nearby companions is currently too small to de- 
termine if there is a correlation between the presence 
of a companion and the hot Jupiter orbital or physical 
properties (or host star type). Furthermore, the cluster 
of hot Jupiters around an orbital period of 2-3 days is 
dominated by transit detections, whereas hot Jupiters 
are quite common on 4-5 day orbits when considering 
only radial velocity discoveries. Further analysis is re- 
quired to quantify the extent to which the hot Juiters 
with close companions are outliers in the period-radius 
space. 

What we know today is that the known hot Jupiter 
nearby companions predominantly orbit interior to the 
hot Jupiter, so the hot Jupiter must orbit at a certain 
distance from the host star so that there is enough room 
in the system for a stable orbit of a smaller planet inte- 
rior to the hot Jupiter. There is the notable exception of 
the WASP-47 system, which also has one exterior com- 
panion. Whether additional exterior transiting compan- 
ions to hot Jupiters exist remains to be seen, especially 
given the fact that all of the currently known close com- 
panions have been discovered via transits, which are bi- 
ased towards planets on shorter orbital periods. With 
the probability of transit dropping with increasing or- 
bital period, the discovery of additional exterior tran- 
siting companions is more challenging. Constraints have 
been placed on the existence of nearby, non-transiting 
planets through the investigation of transit timing varia- 
tions (TTVs) using Kepler data, with only WASP-148 c 
found (Steffen et al. 2012; Hébrard et al. 2020), support- 
ing the idea that exterior companions are indeed perhaps 
intrinsically rarer than interior companions. Analysis of 
comprehensive transit timing searches using TESS data 
(such as that of Ivshina & Winn 2022) could provide 
updated constraints, however. 

Simulations performed by Ogihara et al. (2013, 2014) 
have provided a possible explanation for the dearth 
of external companions. They show that the existence 
of super-Earths exterior to a hot Jupiter would drive 
the hot Jupiter into the host star. Therefore, any hot 
Jupiters with external planets would have been driven 
into the host star prior to their discovery, thus leading 
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Figure 8. The period and radius of all planets in the four hot Jupiter systems with known companion planets (WASP-47, 
Kepler-730, TOI-1130, WASP-132) overlaid on the periods and radii of all planets contained on the NASA Exoplanet Archive. 
The larger, colored markers represent planets in one of these systems, with similar markers corresponding to planets in the same 
system. The markers for WASP-132 b and c are outlined in black. A dashed horizontal line has been included at 8Rg to denote 


the division between hot Jupiters and smaller planets. 


to the “crowding-out of giants by dwarfs”, as Ogihara 
et al. (2013) refers to it. 

In each of the previously known hot Jupiter systems 
with nearby companions, none of the planets are in or- 
bital resonance with each other. The case is the same 
for WASP-132. A small fraction of Kepler multiplanet 
systems are in resonances (Lissauer et al. 2011), whereas 
wide-orbiting giant planets frequently are (Winn & Fab- 
rycky 2015). Expanding the sample of hot Jupiters with 
nearby companions would allow for the comparison of 
the period ratio distribution of hot Jupiters with nearby 
companions to other known exoplanet populations. Sim- 
ilar distributions of period ratios would be expected 
from populations of systems that were assembled in a 
similar manner, and comparing the period ratio distribu- 
tion of hot Jupiters with companions and the occurrence 
of resonances to a population such as super-Earth sys- 
tems could elucidate possible common formation mech- 
anisms. 


7.3. Prospects for Follow-Up 


These hot Jupiters with companions provide a unique 
opportunity to test if spin-orbit misalignment is a con- 
sequence of high-eccentricity migration as theorized. 
Based on the derived parameters of the host star and 
WASP-132 c, we estimate that the semi-amplitude of 
the Rossiter-McLaughlin effect for WASP-132 c is ~0.5 


m s^! (assuming a v sin i ~0.9 km s^! as reported by 


Hellier et al. (2017)). The estimated semi-amplitude of 
the hot Jupiter WASP-132 b is much larger at ^17 m 
s71, Although the prospect of measuring the Rossiter- 
McLaughlin effect for WASP-132 c is quite low, that of 
WASP-132 b can easily be detected with current radial 
velocity precision and may inform our understanding of 
the orbital alignment with the host star's rotation. The 
WASP-132 system is otherwise amenable to additional 
follow-up observations as it is relatively bright at near- 
infrared wavelengths (K, = 9.674). 

Further radial velocity monitoring and a measured 
mass is necessary to confirm the planetary nature of 
this planet. Given that the expected semi-amplitude of 
WASP-132 c is between 3 and 6 ms~!, an instrument 
more sensitive than CORALIE such as HARPS or PFS 
would be ideal to follow up on WASP-132 c to obtain a 
mass measurement and further characterize the system. 
Both of these instruments have 1 ms~! sensitivity or 
better, making a detection of WASP-132 c with either 
instrument possible. 


8. SUMMARY 


In this paper, we present the discovery and validation 
of a companion planet orbiting interior to hot Jupiter 
WASP-132 b. Our investigation and results are summa- 
rized here: 
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1. A ~1.01 d periodic signal with a false alarm 
probability < 1% was detected in the TESS 
photometric data for WASP-132 (TOI-822) with 
both the SPOC pipeline and the Transit Least 
Squares search algorithm. Neither the TESS 
SPOC pipeline nor the vetting software DAVE 
found any immediate false positive indicators. 


2. We refined the system parameters, obtaining a 

host star of R, = 0.753* 002» Ro and M, = 0.782 
+ 0.034 Mo and planets with periods and radii of 
7.13 d and 10.05 + 0.28 Rg for WASP-132 b and 
1.01 d and 1.85 + 0.10 Re for WASP-132 c. 


3. An analysis of archival CORALIE radial velocity 
measurements did not yield a significant detection 
at the 1.01 d period, with a 3c upper limit of 37.35 
Ma on the mass of WASP-132 c. 


4. Using LCOGT ground-based follow-up photom- 
etry, we ruled out NEB signals as the potential 
source of the TESS detection in stars out to 50" 
from the WASP-132 c, and likely ruled out poten- 
tial NEB signals in stars out to 2/5 


5. WASP-132 c is statistically validated as a planet 
with false positive probabilities (FPPs) of 9.02 x 
10-75 and 0.0107 using vespa and TRICERATOPS, 
respectively (see Section 5.2 for further discus- 
sion). 


6. The system is dynamically stable on timescales 
of 100 Myr for planetary and orbital parameters 
within 3c of the best-fit values. 


7. WASP-132 is the second system discovered by 
TESS (and one of only four systems in total) to 
contain a hot Jupiter with a nearby companion 
planet, suggesting that a mechanism other than 
high-eccentricity migration may play a significant 
role in the formation of hot Jupiters. 


'The discovery of WASP-132 c demonstrates the abil- 
ity of TESS to not only find new planets but also en- 
hance our knowledge of those already known. As TESS 
continues its almost-all-sky survey, it will surely reveal 
additional systems similar to WASP-132 which will im- 
prove our understanding of the evolution of hot Jupiter 
systems. With an ever-expanding census of hot Jupiters 
with nearby companion planets, it may even be possible 
to identify sub-populations of these hot, giant planets. 
It is imperative to continue the search for this type of 
system architecture as a larger data set will allow us 
to solidify our understanding of how these rare systems 
form. 
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